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Abstract. We study management of localized modes in two-component (spinor) 
Bose-Einstein condensates embedded in optical lattices by changing interspecies 
interactions. By numerical integration of the coupled Gross-Pitaevskii equations, 
we find three different regimes of the delocalizing transition: i) the partial 
delocalization when the chemical potential of one of the components collapses 
with a gap edge and the respective component transforms into a Bloch state, 
while the other component remains localized; ii) the partial delocalization as 
consequence of instability of one of the components; and iii) the situation where 
the vector soliton reaches limits of the existence domain. It is shown that there 
exists a critical value for interspecies scattering length, below which solutions can 
be manipulated and above which one of the components is irreversibly destroyed. 
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Management of matter waves in optical lattices (OLs) can be implemented in 
different ways, including variation of parameters of the periodic potential in time 
or change of the nonlinearity with help of the Feshbach resonance. In particular, 
by changing the amplitude of the lattice one can achieve the delocalizing transition, 
which consists in transformation of an initially localized mode, also referred to as a 
gap soliton, into a spreading out atomic wave packet. Such a transition was predicted 
for two- and three-dimensional single-component Bose-Einstein condensates (BECs) 
[TJ H] and can be observed by adiabatically decreasing the amplitude of the periodic 
potential below some value, with subsequent restoration of the initial lattice profile. 
If the amplitude of the lattice reaches some critical value, the initially localized wave 
packet cannot be recovered at the end of the experiment and one observes a spreading 
out condensate. In a recent paper [3] , it has been shown that the delocalizing transition 
can be observed also in one-dimensional (ID) systems, with help of interplay between 
linear and nonlinear lattices. This suggestion was based on the fact that the transition 
occurs in a system where small amplitude gap solitons cannot exist [5], i.e. where a 
minimal atomic number is necessary for the creation of a localized mode. From the 
physical point of view, the transition is related to modulational stability of the Bloch 
states bordering either the semi-infinite gap in the spectrum of Bogoliubov phonons 
or a gap induced by the periodic potential [3j. 

The described results naturally rise questions: Whether other ID systems allowing 
for delocalizing transition exist? and Whether other physical mechanisms resulting 
in delocalizing transition exist? A simple positive answer to the first question readily 
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follows from the fact that the quintic nonlinear Schrodinger equation with a periodic 
potential requires a critical number of particles for the creation of a gap soliton [1], 
and thus can undergo the delocalizing transition. The answer to the second question 
is also positive, and moreover it offers another ID system allowing for a new type 
of delocalizing transition. To explain this, we recall that new factors affecting the 
existence of gap solitons appear when one considers multicomponent systems in 
periodic potentials, and, in particular, binary mixtures or spinor condensates loaded in 
OLs: a necessary condition for existence of the modes is that the chemical potentials 
of both components must belong to the respective phonon gaps [5]. This condition 
may be violated being subjected to change of the system parameters, thus leading to 
destruction of two-component gap solitons. Generally speaking, in such a situation 
one does not necessarily expect complete delocalization of the condensate, because 
one of the chemical potentials may still belong to a gap, while the other one reaches 
the gap edge. Instead, in this situation, one may observe a partial delocalizing 
transition, where one of the components spreads out, while the other one stays 
localized. Detailed description of this phenomenon constitutes the main goal of the 
present communication. 

Before going into details, we mention other aspects of the formulated problem. 
First, transformation between two different states of a condensate represents not only 
fundamental but also practical interest, as it allows for selective and thus controllable 
separation of the initially localized components of mixtures of BECs. Its mathematical 
relevance is determined by the fact that the transition allows for experimental study 
of the limits of the existence of localized states. 

Second, on the basis of preceding studies, the most direct approach to obtain the 
delocalizing transition would be the change of the parameters of the OL. However 
(this was checked by our extensive numerical simulations), the response of a binary 
mixture to such change is very similar to the one displayed by one-component BECs, 
and does not lead to the delocalizing transition. In particular, decreasing (increasing) 
of the potential results in expansion (narrowing) of vector gap solitons (c.f. [BE]). 
Essentially new phenomena can be found when management is performed by means 
of changing the interspecies interactions (for one-component BECs, the delocalizing 
transition induced by nonlinear management was suggested in |3J). 

Third, it was also verified numerically that the delocalizing transition reveals 
qualitatively the same behavior for binary mixtures of two different species of atoms 
and for spinor condensates. Therefore, in the present work we concentrate on the 
more simple dynamics of a spinor condensate induced by change of the interspecies 
interaction. 

Finally, since adiabatic change of the parameters in practice represents a 
particular type of management of mixtures of BECs, it is relevant to mention that, 
in Ref. [6], it was shown that interplay between the OL potential and the scattering 
length can assist building two-component gap solitons. 

The model and numerical method. A diluted binary mixture of = 1/2 spinor BECs 
in the mean-field approximation in a quasi-lD geometry is described by the coupled 
Gross-Pitaevskii (CP) equations 



y cos(2a;)V'j + {xjl^jf + xl^s^jf) ^j- 



(1) 



dt 
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Here, j = 1,2, and scaling of the variables is the same as in |5j. In particular, 
in the chosen units, the effective number of particles of component j is calculated 
as Nj = J \ijjj\'^dx and is related to the real number of atoms, Afj, by the formula 
Nj = gMj, where g ~ 10"'^ 10^'* in the real experimental situation. The amplitude 
of the OL, V, is measured in the units of the recoil energy. The coefficients Xj and x 
characterize two-body interactions. 

Since we are interested in a situation where x = x(0 adiabatically changing 
in time, stationary localized solutions of Eq. ([1]) play a prominent role. We thus 
make an ansatz il}j{x,t) = ^'j(x)e~*'^J*, where fij is the chemical potential of the j*'' 
component, and substituting in ([T]) we obtain {j — 1, 2) 

" - l^cos(2x)*, + (x,|*,f +xl*3-,f ) (2) 

A diversity of localized modes supported by Eq. (U) possessing a given symmetry 
were described in [5] . Here we focus only on the particular case of symmetric (even) 
modes. Moreover, to simplify the construction of the initial states (we do this by 
the shooting method), as well as for their systematic exploration, we restrict the 
consideration to the case where (xo — Xi)(Xo — X2) > 0, where xo is defined as 
Xo = x(i — 0)- In. this case, the simplest initial profile of a localized mode is 
characterized by equal chemical potentials /^i = /X2 = /xq. Then = a~^$, where 
aj = VlXo ~ X1X2I/IX0 — X3-jl a-nd the real function $ solves the equation 

^$ = - V cos(2a;)$ + cr$3 (3) 

with a =:sgn[(x§ — XiX2)(xo ~ X3-i)] (due to the above constrain, it does not depend 
on j). A relation among the numbers of particles is given by / ~ N = aiNi = 

alN2. 

For each of the components, the OL results in a band spectrum, and chemical 
potentials must belong to the respective gaps. We further restrict our analysis to /xq 
belonging either to the semi-infinite or to the first lowest gap. To be specific, we will 
also fix y = 1, which does not affect qualitatively the phenomenon. 

As we investigate the behavior of localized modes subjected to the adiabatic 
variation of Xi we first employ a Newton-Raphson (NR) method for constructing the 
branches of solutions, which "depart" from the given ^E*-,- and are constrained to fixed 
values of Ni^2- When the strength of the interspecies interactions is changed, the 
chemical potentials of both components arc changed as well, i.e. /ij = ^j(x)i provided 
that /ii(xo) = M2(xo) = Mo, and can be computed using the formula (j = 1, 2) 



9ti ^ 



dx 



V cos{2x)\'ipj 



|2 



+ xM,r + x\^3-,Vm')dx. (4) 

This formula is also used to calculate "instant" chemical potentials in the temporal 
adiabatic dynamics, when x = x(0' which is studied using the procedure as follows. 
For a given xo, we start with the solutions — aj'^f^ obtained from Eq. ([3]). For the 
sake of definiteness, variation of the interspecies interaction is chosen as 

X{t) = (XO - Xm) COS^ ^ 

with Xm = x(^/2) and T being the total time of simulations. In other words, we first 
increase (decrease) the nonlinearity until the maximal (minimal) value Xm and then 
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decrease (increase) it until the initial value xo- Behavior of this system is investigated 
by a direct integration of Eq. ([1]), with adiabatic variation of as in Eq. 

For management to be possible, the initial solutions must be experimentally 
feasible. Therefore we test the dynamical stability of the initial 5*1,2, by perturbing 
it with some addenda (in the presented results it is given by 0.05 cos(6.5(a; + 0j))\E'j, 
where 0i = 0.0 and 02 = 0.3) and numerically integrating Eq. ([1]), with constant 
X = XO- 

Repulsive two-body interactions. First we consider a binary mixture with all repulsive 
interactions. To be specific we choose xi = 1-Oj X2 = 0.5 and xo = 2.0, what 
corresponds to ai = 1.53, a2 — 1.87 and a = \. We start with the existence of 
the stationary fundamental vector soliton for different using the NR method, as 
explained above. The chemical potentials /ii.2(x) ^"^^ depicted in Fig. [T] (the left top 
panel), where one can observe two important features. 




Figure 1. The chemical potentials vs nonlinear coefficient x (l^ft top panel), for 
A'^l ~ 1.951 and N2 ~ 1.303. Dotted lines correspond to the limits of the first 
energy gap (fi G [0.471; 1.467]). fii* is the chemical potential of the gap soliton 
of component 1, with A^i 1.951, computed at 'I'2 = (see also Fig. |2)l. In 
the panels A, B and C, the profiles of stationary modes SPi (solid lines) and ^2 
(dashed lines) in the points A, B and C are shown. The initial chemical potential 
is fio = 1.3. 

First, even infinitesimal change of x = Xo + ^X results in splitting of the chemical 
potentials: fij — fj^o -\- Sfij with 7^ (5/^2. Second, the chemical potential of the 
second component monotonously grows with x approaching the top of the first energy 
gap. The latter suggests that, at some x — Xcr between points B and C (see the left 
top panel of Fig. [1]), ^2 collapases with the gap edge. For x > Xcr, the vector soliton 
with given numbers of particles in each component does not exist and thus, according 
to the arguments exposed above, Xcr must determine the critical value above which 
the delocalizing transition occurs. Moreover, one can predict that this is a partial 
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delocalization because the chemical potential of the first component still belongs to 
the gap and thus the respective component persists. 

Although we use a rather large spatial domain of calculation x g [— 487r, 487r], the 
provided analysis of the stationary solutions is not yet conclusive about the existence of 
the modes, because of boundaries possibly affecting the expanding second component. 
Indeed, from the profiles of shown in Fig. [1] one observes tendency of the second 
component to the Bloch state of the upper band edge as x exceeds Xcr- To check 
numerically whether the delocalization is authentic, i.e. that the second component 
does not represent, say, a small amplitude gap soliton, affected by the boundaries, 
we recall that |5'2p becomes small enough when x tends to Xcr- Hence the first 
component approximately solves Eq. ([2]) with = 0. That is why, in Fig. [2l we 
constructed the lowest branch of the localized modes of the first component in the 
absence of the second component. For a given number of particles Ni, we found the 
chemical potential /i* and constructed the one-component localized mode. 




Figure 2. (a) The number of particles of the first component along the first 
energy gap {fj, G [0.471; 1.467]), with '3/2 = 0. Horizontal thin line represents 
A^i 1.951 (i.e. the number of particles of the first component in the mixture 
analized in Fig. [T]l. (b) The solid line represents the profile of ^"1 from Fig. [TJi;. 
The dashed line is obtained as a solution of Eqs.l[2]l, with ^2 = and A'l = 1.951 
(the lines are indistinguishable on the figure scale), (c) The solid line corresponds 
to 4*2 taken at the point C of Fig. [T] compared with the Bloch wave (the dashed 
line) calculated at the lower edge of the second band {/i = 1.467) and having an 
amplitude matching the solution >I'2- 

In Fig. [TJleft top panel), one can see that, indeed, fii — > /ii* (see the dashed 
line) in the domain where /i2 collapses with the second band. Direct comparison 
of the explicit shapes of the localized modes in the single-component BEC and the 
first component in the binary mixture is shown in Fig. [IJb), illustrating that the 
distribution of the first component is transformed into the one-component localized 
mode. In the same point, the second component becomes an extended Bloch state, as 
shown in Fig. [2jc) . 

Finally, we turn to the dynamical experiment. The results are shown in Fig. [3l 
Scanning values of Xm between those corresponding to the points B and C in Fig. [U 
one observes complete restoring of the vector soliton if x does not reach its critical 
value, which is found to be Xcr ~ 2.6, while partial delocalization of the binary mode 
is observed when Xm > Xcr (as predicted by the stationary analysis). 

Additionally, we explored the case with all attractive interactions {x,Xj < 0)- 
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Figure 3. In (a), time dependence of the intraspecies interaction is shown. Solid 
and dashes lines correspond to Xm = 2.5 and Xm = 3. Dynamics of the densities 
of the first component in (b), (d) and of the second component in (c), (e) are 
presented. Xm = 2.5 for [(b), (c)], and Xm = 3.0 for [(d), (e)]. Here, = 1-3 and 
XO = 2.0. 



Symmetric solutions are now located in the semi-infinite gap and, as x is adiabatically 
increased, the chemical potentials of both components grow, approaching the gap edge. 
However, we did not observe the delocalizing transition for negative values of x- It 
occurred only for x > 0, i.e. when Xj ■ X becomes negative. Therefore, next we focus 
on this case. 



Cases ■ X < 0- Now we consider a binary mixture when intra- and interspecies 
interactions have different signs. To be specific, we investigate the case where 
Xi = —1.0, X2 — —0.5 and x 0- In such a case, both chemical potentials 
belong to the semi-infinite gap. Starting with uncoupled components, xo = 0.0, 
what corresponds to ai = 1.0, a2 = 0.71 and a = —1, we adiabatically increase x- 
Behavior of the chemical potentials fJ.i^2ix) i^i the stationary problem and snapshots of 
the profiles of the components of vector solitons at different values of x are presented 
in Fig. m 

As in the previous cases, change of x results in splitting of chemical potentials 
Hj. Also, there exists a maximum strength for the repulsive interspecies interaction, 
X — Xex, such that the fundamental vector soliton of the stationary problem can be 
found for x < Xex- This point corresponds to the bifurcation point between two vector 
solitons whose profiles are represented in Fig.|4]B and C. Thus, contrary to the case of 
all repulsive interactions, now Xex is a bifurcation point where /ii,2 belong to the gap 
but where 8111^2/ dx = 00. For the numbers of atoms as in Fig.[?l it was numerically 
found that Xex ~ 0.0975. Passing to the dynamics, we integrated the GP equations 
([1]). For different values of Xm below the bifurcation point, Xm < Xex, we found that 
the vector soliton can be restored to its initial shape, as shown in[5[(b),(c)]. However, 
if X overcomes the critical value Xex (see Fig. (S^d), (e)]), one component (here ^'i) 
breaks into two repelling wave packets moving outwards the center, while the other 
component (here ^'2) transforms into a stable single-component localized mode. The 
observed behavior can also be interpreted as a kind of delocalizing transition. Its 
physical origin and, naturally, its manifestation are different from the ones reported 
in Fig. H 

Additionally, instability of the vector soliton was observed for several values of 
Xm < Xex ■ AH such valucs are located in a region of the branch fj,2 (x) between its 
inflection point and Xex- We attribute this behavior to the existence of domains of 
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Figure 4. The chemical potentials vs nonlinear coefficient x (left top panel). The 
dotted line corresponds to the upper edge of the semi-infinite gap. In A, B and C, 
the profiles of stationary modes 'I'l (solid lines) and \E'2 (dashed lines) are shown. 
Here, A^i « 0.996 and N2 ~ 1.993. The initial chemical potential is /iq = —0.25. 
When X approaches Xex, solutions B and C bifurcate. 



instability along the branch /i2(x)- From Eq. ([5]), we have that the velocity of x(i) at 
the point Xm is zero (xiXm) = 0), therefore a mode spends long enough time in the 
vicinity of Xm for the instability to develop. This is what happens in Fig. [5[(f), (g)], 
where one of the components of the vector soliton (here ^'i) looses its stability, breaks 
from the second component and starts to move away from the center, while the second 
component maintains its shape. In the case where Xm > Xex, presented in Fig. [5j(d), 
(e)] , although we pass in the instability region, we do not observe instability. This can 
be explained by that fact that the time spent by the vector soliton in the instability 
region was not enough for instability to start developing. 

Conclusions We have studied the management of coupled gap solitons of a two- 
component spinor BEC embedded in an OL, by means of the adiabatic variation 
of the interspecies interaction. When all nonlinear interactions are repulsive, one 
of the components can be delocalized if subjected to a sufficiently large change of 
the interspecies scattering length, as its chemical potential collapses with the energy 
band edge, while the other component remains localized. On the other hand, if 
the intra- and inter-species interactions have different signs, partial delocalization 
is achieved when solutions reach limits of the existence domain. In this case, the 
manifestation of the phenomenon is different, as one of the components breaks in two 
gap solitons propagating outwards each other. Finally, we observed one more cause 
of the delocalizing transition: the instability of a vector soliton which develops more 
rapidly than the change of the scattering length occurs. 

H.A.C. acknowledges support of the FCT through the grant SFRH/BD/23283/2005. 
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Figure 5. In (a), the time dependence of the inter-species interaction is shown. 
Sohd, dash-dotted and dashed lines correspond to Xm = 0.02, Xm = 0.08 and 
Xm = 0.15. Dynamics of the densities of the first component in (b), (d), (f) 
and of the second component in (c), (e), (g) are presented. Xm = 0.02 for 
[(b), (c)], Xm = 0.15 for [(d), (e)], and Xm = 0.08 for [(f), (g)]. Initial parameters 
are /io = —0.25 and xo = 0.0. 
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